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Recently, it has been reported that a candidate for a quark star may have been observed. In this 
article, we pay attention to quark stars with radiation radii in the reported range. We calculate 
nonradial oscillations of /-, w- and u>n-modes. Then, we find that the dependence of the /-mode 
quasi-normal frequency on the bag constant and stellar radiation radius is very strong and different 
from that of the lowest wn-mode quasi-normal frequency. Furthermore we deduce a new empirical 
formula between the /-mode frequency of gravitational waves and the parameter of the equation 
of state for quark stars. The observation of gravitational waves both of the /-mode and of the 
lowest u>ii-mode would provide a powerful probe for the equation of state of quark matter and the 
properties of quark stars. 
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I. INTRODUCTION 



Since the 1980s, there have been many studies about objects known as quark stars. If the true ground state of 
matter is not the nuclear state but a quark state as Witten conjectured there may exist objects composed of 
quark matter. For simplicity, the bag model equation of state (EOS) has often been used for studying quark matter. 
' Lattimer and Prakash calculated stellar properties with this EOS Q- In reality, it is quite difficult to determine the 
EOS of quark matter by experiments on the ground, because quark matter will appear at extremely high density. In 
this respect, the observation of astrophysical compact objects may be a unique means to determine the equation of 
state for cold high-density matter 

Recently, the radiation radius of the compact star RXJ185635-3754 was estimated to be very small (3.8km < i?^ < 
8.2km) based on X-ray observations by Chandra The radiation radius is related as Rqo = R/ \/l — 2M/ R to the 
areal radius R and mass M of the star. This compact star is too small to construct with the normal EOS's which are 
| currently adopted for neutron stars Q. This star should be a candidate for a quark star. Inspired by this observation, 
O" 1 Nakamura proposed several scenarios about the formation of quark stars @ . 

However, it should be noted that there exist many subsequent articles which question the whole observation. For 
example, Pons et al. pointed out, by employing the data of not only X-ray observation but also UV/optical observation, 
that if RXJ185635-3754 is represented by the simplest uniform-temperature heavy-element atmospheric model, this 
• i-H , compact star has Roo « 8km, M 0.9M Q and R m 6km 7J. Though this radiation radius is included in the range 
^ suggested by [Bj], these properties can not be constructed, even if the quark matter EOS is used. 

If a compact object oscillates for some reason, gravitational waves are emitted from the object. The trigger of 
the oscillation may be a nonsphcrical supernova explosion, coalescence with another compact object, violent mass 
influx due to possible instability of the surrounding accretion disk, and so on. The oscillations damp out because the 
gravitational waves carry away oscillational energy. So these oscillations are called quasi- normal modes (QNMs). If 
the gravitational waves from a compact object are detected, we can get information about the source object. 

Gravitational wave observation projects, such as LIGO @, VIRGO 9], GEO600 [13] and TAMA300 [H], are 
making remarkable progress in sensitivity. Among them, TAMA300, GEO600 and LIGO are now in operation. If 
highly sensitive gravitational wave detectors are available, we may obtain a large amount of data for frequencies and 
damping times of gravitational waves emitted from compact objects. The systematic study of gravitational wave modes 
leads to the so called "gravitational wave asteroseismology" , which has been initiated and being presented by various 
authors (see 0, Hi Hi El ) • If we have a good empirical formula for gravitational wave oscillational frequencies 
and damping times as functions of stellar properties, it will be very useful in obtaining information for source stars 
from gravitational wave observations. In this context, Andersson and Kokkotas calculated QNMs with various EOSs 
and proposed an empirical formula between the QNMs and properties of neutron stars 0>EI- Kokkotas, Apostolatos 
and Andersson improved this empirical formula by one which included the relevant statistical errors . These works 
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treated the polar modes of oscillation and the results are well summarized in [Til HtH . Benhar, Berti and Ferrari 
calculated axial modes and showed that the axial mode gives more direct and explicit information on the stiffness of 
EOS compared to the polar mode All these works are for normal neutron stars, i.e., assuming normal equations 
of state for nuclear matter. 

On the other hand, Yip, Chu and Leung studied nonradial stellar oscillations of quark stars whose radii are around 
10 km 0). However, the reported radius of the star RXJ185635-3754 is much smaller than the radii of the stellar 
models examined by Yip, Chu and Leung. In this paper, we calculate nonradial oscillations of quark stars whose 
radiation radii are in the range of 3.8km < i?^ < 8.2km, and examine the relation between the QNMs of quark stars 
and the EOS of quark matter. Our main concern in this article is whether one could distinguish the EOS of quark 
matter by the direct detection of the gravitational wave. 

QNMs are classified into two classes. One is a class of fluid modes, which mainly couple with the stellar fluid, 
while the other is one of spacetime modes which are connected with spacetime oscillation, which couple mainly to 
metric perturbations. Non-rotating stars admit fluid modes only for polar perturbations. The most widely studied 
fluid modes are /-, p- and (/-modes 0,0]. The /-mode is the fundamental mode. There is only one /-mode for 
each I. The p-mode is the pressure mode which arises from fluid pressure. The g-mode is the gravity mode which 
is caused by buoyancy due to density discontinuity and/or the temperature gradient of stars and other factors. The 
damping rates of these fluid modes are very small compared with the frequency. The spacetime modes include w- 
and wit- modes [53, ■ The damping rate of w-modes is comparable to the frequency. For each stellar model, there 
are found one or a few u>n-modcs, whose complex frequency is located near the imaginary axis. Here we deal with 
only polar perturbation because our main concern is the relationship between the gravitational waves and EOS. For 
simplicity, we pay attention here only to /-, w- and wn-modes. Recently, Kojima and Sakata calculated the /-mode 
quasi-normal frequency for quark star models, and pointed out the possibility of distinguishing between quark stars 
and neutron stars by detecting both the /-mode frequency and damping rate p^ . 

The plan of this paper is as follows. We present basic equations to construct quark stars and show the properties 
of quark star structure in Sec. [n] I n Sec. 11111 we present the equations and method to determine the QNMs for 
spherically symmetric stars, show numerical results for the quasi-normal frequencies of quark-star models constructed 
in Sec. [HI and discuss their implications. Moreover, we apply the empirical formula proposed by Kokkkotas et al. 
to quark star models and deduce a new empirical formula from our numerical results between the /-mode frequency 
of the gravitational wave and the parameter of EOS in Sec. II VI We conclude in Sec. [V] We adopt the unit of 
c = G = 1, where c and G denote the speed of light and gravitational constant, respectively, and the metric signature 
of (-,+,+,+). 

II. STRUCTURE OF QUARK STARS 

We consider a static and spherical star. For this case, the metric is described by 

ds 2 = -e 2<s> dt 2 + e 2A dr 2 + r 2 (dO 2 + sin 2 9d(f> 2 ) , (2.1) 
where <£>, A are metric functions of r. The mass function m(r) is defined as 

m(r) = -r (l - e~ 2A ) . (2.2) 

This mass function satisfies 

dm n 

— = ^r 2 p, (2.3) 

where p is the energy density. The equilibrium stellar model is derived by solving the Tolman-Oppenheimer-Volkoff 
equation, 

dP _ (p + P) (to + 47rr 3 P) 
dr r (r — 2m) 

where P is the pressure of the fluid, and the potential 4> is given by 

d<f> (to + 47rr 3 P) 
dr r (r — 2m) 

For solving these equations, we need an additional equation, i.e., the equation of state. 



(2.5) 
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Though a quark star may be either a light quark star made up of u and d quarks or a strange quark star composed 
of a mixture of u, d and s quarks, we construct simple models by using the bag model EOS, which neglects the masses 
of it and d quarks. This bag model EOS is characterized by three parameters: the strong interaction coupling constant 
a c , the bag constant B and the mass m « of the s quark. The dependence of stellar properties on the bag constant is 
much stronger than that on a c and m s [23J. In this paper we use the bag model EOS derived for a massless strange 
quark, 

P=i(p-4B). (2.6) 

The bag constant B is a positive energy density, which corresponds to latent heat. The stellar radius R is defined 
as the position where the pressure is zero. At the surface, the density profile is discontinuous. Although the value of 
the bag constant, which is a phcnomcnological parameter, should be determined by the underlying strong interaction 
dynamics, it is difficult to determine this value from our present understanding of QCD. However, if strange quark 
matter is the true ground state, the energy of this true ground state per particle would not be over the nucleon mass 
939 MeV at baryon density for zero pressure matter. This constraint implies a maximum value of the bag constant 
B as follows 01: 

B < S max = 94.92 ( 1 - — ) MeV fm~ 3 . (2.7) 



For definiteness, we use three values of the bag constant, B = 28.9, 56.0 and 94.92 MeV fm~ 3 . Although bI( 4 = 145 
MeV or B st — 57.8 MeV fm -3 has been often used for the study of quark-gluon plasma so far, we adopt B = 56.0 
MeV fm~ 3 to compare with the results in [l^. The other values of the bag constant are the maximum value for 
a c — and half the standard value B st . We show the mass M of quark stars as a function of the central density p c 
for these three values of the bag constant in Fig. Qfa), the relations between M and R and between M and R^ in 
Fig. IHb), and the relation between "average density" p = 3M/4ttR 3 and R^ in Fig. [He). These figures are plotted 
for 4B < p c < 5.0 x 10 15 g/cm~ 3 . The radiation radius may be determined from X-ray observations such as p|. We 
pay attention to the stars whose radiation radii are in the range of 3.8 — 8.2 km. In this range of radiation radius, we 
pick up three values, R^ = 3.8, 6.0 and 8.2 km. 

On the other hand, Nakamura argued that the mass of the compact star reported in p| should be roughly 0.7M Q 
in order to account for its observed X-ray luminosity If we use a value smaller than -B max , however, the mass of a 
quark star whose radiation radius is in the range of 3.8 — 8.2 km is far below 0.7M Q as seen in Fig. db). Therefore 
he pointed out the possibility of adopting the value B = 471.3 MeV fm -3 , which is rather greater than -B max - We 
consider a quark-star mass of M = 0.7M Q for this value of the bag constant to compare with other models. Also for 
the case B = 471.3 MeV fm -3 , we plot the relation between M and p c in Fig. d a ); the relations between M and R 
and between M and Roo in Fig. [Hb), and the relation between p and i?oo in Fig. [He). In this case, each plot is for 
4B < p c < 5.0 x 10 16 g/cm~ 3 . The properties of our quark-star models are tabulated in Table [U 



III. NONRADIAL OSCILLATIONS OF QUARK STARS 

A. Method 



The QNMs are determined by solving the perturbation equations with appropriate boundary conditions. The metric 
perturbation is given by 



a - a {B) 



(3.1) 



where gffl is the background metric of a spherically symmetric star l|2.1|l . We have applied a formalism developed by 
Lindblom and Detweiler |25| for relativistic nonradial stellar oscillations. In this formalism, h^y for polar perturbations 
is described as 
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where H,H\, and K are perturbed metric functions with respect to r, and the components of the Lagrangian dis- 
placement of fluid perturbations are expanded as 
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where W and F are functions of r. 

Then the perturbation equations derived from Einstein equations are given by 
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(3.7) 
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(3.10) 
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In deriving the above equations for perturbations, we assume a perfect fluid. Rigorously speaking, it will not be valid 
because of the existence of matter viscosity. However, little is known about the viscosity of cold quark matter. Here 
we omit the viscosity as a possible approximation. In Eq. (|3.8f) . 7 is the adiabatic index of the unperturbed stellar 
model, which is calculated as 
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(3.12) 



where the bag model EOS (|2.6[1 is used in the second equality. The set of Eqs. I|3.6|) — (I3.9|l is a set of differential 
equations connecting the variables H\, K, W and X, and Eqs. (|3.1l)|l and (|3.11|l are the algebraic equations for the 
variables H and V. The perturbation equations outside the star are described by the Zerilli equations. By imposing 
boundary conditions such that perturbative variables are regular at the center of the star, the Lagrangian perturbation 
of pressure vanishes at the stellar surface, and the gravitational wave is only an outgoing one at infinity, one can reduce 
this to an eigenvalue problem. The boundary condition at the stellar surface is X ~ 0, because X = —r _ 'e*AP, 
where AP is Lagrangian perturbation of pressure. Furthermore, we set the term e~®X /7P in Eq. I|3.8|l to zero at the 
stellar surface, because 7P = 4P /3 at r = R. For the treatment of the boundary condition at infinity, we adopt the 
method of continued- fraction expansion proposed by Leaver [26| . The details of the determination of quasi- normal 
frequencies are given in |l9j . 
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B. Numerical Results 

1. f-mode 

We plot the complex frequencies of the /-mode for I = 2 in Fig. In this figure, the squares, triangles and circles 
denote the cases B = 28.9 MeV fm" 3 , B = 56.0 MeV fm~ 3 and B = 94.92 MeV fm^ 3 , respectively. In each set, the 
upper, middle and lower marks correspond to the stellar models of Roc = 8.2 km, i?^ = 6.0 km and i?oo = 3.8 km, 
respectively. The labels in this figure correspond to those of the stellar models in Table [I] 

As can be seen in Fig. EI a), the /-mode frequency depends strongly on the bag constant but very weakly on the 
stellar radiation radius. This result agrees with that of Kojima and Sakata j2^. Therefore, if the radiation radius of 
a quark star is determined by observation, even though the value of the radiation radius is somewhat uncertain, we 
can directly obtain the value of the bag constant by detecting /-mode QNMs (see Sec. IIV|I . Furthermore, we can 
restrict the mass of the source star by employing the observational value of the radiation radius. 

The reason why the /-mode frequency depends strongly on the bag constant is understood as follows. The /-mode 
frequency depends strongly on average density but very weakly on the EOS . In the case of quark stars with small 
radii, the average density is almost determined by the bag constant alone, as can be seen in Fig. (He). Thus, the 
/-mode frequency depends strongly on the bag constant. 

Moreover, we can see that the damping rate of the /-mode depends on the stellar radiation radius for each bag 
constant, which also agrees with the result by Kojima and Sakata [23. So, if we get data on the damping rate of the 
/-mode by observation, we can determine the stellar radiation radius independently of X-ray observations. However, 
one must keep in mind that the measurement of frequency is far more accurate than that of damping rate, because 
the estimated relative error in frequency is about three orders of magnitude smaller than that in damping rate [l4| . 

The /-mode frequency of a quark star with M = 0.7M© and B = 471.3 MeV fm -3 (ATI) is plotted in Fig. Efb). In 
this figure, results for the case B < i? max are also plotted for comparison. Compared with the /-mode frequencies for 
B < B max , the plot for B = 471.3 MeV fm -3 is relatively far away in phase space. This is mainly due to the large 
difference in the average density of quark stars. 

2. w- and wu-modes 

The calculated complex frequencies of w- and u>n-modes for I = 2 are plotted in Fig. In these figures, the filled 
and non-filled marks correspond to w- and tun -modes, respectively. Figures E{a), (b) and (c) correspond to B = 28.9, 
56.0 and 94.92 MeV fm -3 , respectively. In each figure, the upper, middle and lower sequences correspond to the stellar 
models of = 3.8 km (Al, Bl, CI), R^ = 6.0 km (A2, B2, C2) and R^ = 8.2 km (A3, B3, C3), respectively. 
In Fig. Eld), we plot the complex frequencies of w- and wu-modes of the star of O.7M with B — 471.3 MeV fm -3 
(VI). 

In order to obtain the tendency of w- and wu-modes with changing compactness, we calculate the w- and u>n-modes 
of the stellar models whose compactness is much larger than that of the stellar models depicted in Figs. Ofa), (b) 
and (c). From this calculation, we find two wu-modes for the stellar model of higher compactness, for example, 
for p c — 1.4 x 10 15 g/cm -3 with B = 56.0 MeV fm -3 . In this case, the complex frequencies of w- and wu-modes 
appear similarly as in Fig. Old). Furthermore, we see that as the compactness of quark star gets smaller, the complex 
frequencies of the wn-modes shift toward the imaginary axis, while both the frequency and damping rate of w-modes 
become slightly larger. If the compactness is very small, the frequencies of «;n-modes can be too small to find except 
for the lowest wn-mode which has the largest frequency among all wn-modes. That is why Figs. Efa), (b) and (c) 
are quite different from the results obtained by Yip, Chu and Leung [||| in the location of wu-modes. 

We pick up the lowest wn-mode and plot in Fig. 0] for each stellar model. The frequency of the lowest wn-mode 
depends strongly on the bag constant. In this figure, the labels correspond to those of the stellar model in Tabled 
Thus, we can get information about the bag constant or the stellar radiation radius by observing the lowest iun-mode 
frequency. Moreover, we could obtain a more stringent restriction by using the damping rate, although it seems more 
difficult to determine the damping rate by observation as Andersson and Kokkotas mentioned |12| . 

On the other hand, in Fig. |3| we see that both the frequency and damping rate of w-modes have little dependence 
on the bag constant. Therefore it is difficult to determine the value of the bag constant by employing the observation 
of w-modes. In fact, the w-mode frequencies are all far beyond the frequency band for which gravitational wave 
interferometers such as LIGO, VIRGO, GEO600 and TAMA300 have good sensitivity $2%. 

As seen in Fig. Old), because the compactness of model VI is much larger than that of the stellar models constructed 
with B < -Bmax, the appearance of the complex frequencies of the w- and wn-modes is quite different from that for 
the case of B < B max . In Fig^J we also plot the wn-mode complex frequency closest to the imaginary axis among 
the wu-modes previously found for model VI. From this result, we see that if the bag constant is too large, or is as 
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large as the one adopted by Nakamura, it may be difficult to distinguish between the lowest wn-mode for the case 
of B < B max and one which is closest to the imaginary axis for the case of B — 471.3 MeV fm~ 3 . In respect of 
this problem, even if we cannot distinguish between the lowest wu-mode for the case of B < B max and one which is 
closest to the imaginary axis for the case of B — 471.3 MeV fm~ 3 , it may be possible to discriminate between these 
two kinds of wn-modes by employing the observed /-mode frequency, because the plot of /-mode for B = 471.3 MeV 
fm~ 3 is relatively far away in phase space, compared with the /-mode frequencies for B < B max - 



IV. EMPIRICAL FORMULA 



Kokkotas, Apostolatos and Andersson constructed the empirical formulas for /-mode of neutron stars [T^; 
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We apply these empirical formulas to quark star models, and show the results in Table [H] As seen these results, 
though these formulas are deduced by employing various EOS including the realistic one for neutron stars, it is found 
that a simple extrapolation to quark stars is not very successful. We try to construct the alternative formula for 
quark star models by using our numerical results for /-mode QNMs, but we can not do this using the same method 
as Kokkotas et al., because some special relation between /-mode QNMs and the stellar properties can not be found. 

In the previous section, we see that the /-mode frequency depends strongly on not the stellar radiation radius but 
the bag constant. So, in order to get the relationship between /-mode frequency Re(ui) and the bag constant B, we 
calculate /-mode QNMs emitted by six more stellar models in addition to the above nine stellar models (A1^A3, 
B1~B3, C1~C3). The radiation radii of these six added models are fixed at R^ = 3.8, 6.0 and 8.2 km, and the 
adopted bag constants are B — 42.0 and 75.0 MeV fm -3 . The results of calculation are plotted in Fig. and listed 
in table ITTll In Fig. squares, triangles and circles correspond to the stellar models whose radiation radii are fixed 
at 3.8, 6.0 and 8.2 km, respectively. From our numerical results, we find the following new empirical formula between 
the /-mode frequency Re{uj) and the bag constant B: 



Rejuj) 
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0.5252 



(4.3) 



This empirical formula is also shown with a dotted line in Fig. |SJ In Table ITTT1 we listed the bag constant which is 
obtained by substituting the given /-mode frequencies into the empirical formula (|4.3() . and its deviation from the 
true value. Here we calculate the deviation as 



B-B e 
B 



(4.4) 



where B and B emp are the bag constant we adopted for calculation of QNMs and the one given by substituting the 
/-mode frequency into the above empirical formula, respectively. We see that this empirical formula is very useful, 
because it is possible to determine the bag constant precisely by employing the observation of /-mode frequency, even 
if the star has some range of radiation radii. 



V. CONCLUSION 



We calculate the nonradial oscillations of quark stars whose radiation radii are in the range of 3.8 < Roc 8.2 km, 
paying particular attention to the /-, w- and wn-modes of quark stars constructed with the bag model EOS. 

We adopt four values for the bag constant and construct ten stellar models. We find that because the frequency 
of the /-mode strongly depends on the bag constant, we can restrict the bag constant by observations of the /-mode 
frequency. Moreover, since the damping rate of the /-mode depends on the stellar radius for each bag constant, we 
can restrict the stellar radius by detailed observations of the /-mode damping rate. 

On the other hand, both the frequency and damping rate of the lowest wu-mode also strongly depend on the bag 
constant and the stellar radius. Therefore we may be able to get information about the properties of quark stars or 
the EOS governing quark matter from observations of the lowest uin-mode frequency. If the damping rate of this 
mode were obtained, we could derive more stringent constraints on the properties and EOS parameter. 
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Furthermore, we deduce a useful empirical formula between the frequency of /-mode and the bag constant. By 
using this relation, if the /-mode frequency is detected, we can determine the bag constant even if the radiation radius 
of the source star is not determined precisely. Then we can determine the stellar radiation radius by employing the 
observation of the /-mode damping rate. 

The dependence of the /-mode frequency on the bag constant and radiation radius is different from that of the 
lowest wn-mode frequency. Therefore, the lowest u>n-mode QNM can help us to decide the bag constant and/or the 
stellar radiation radius. 

We also calculate the /-, w- and wn-mode complex frequencies for a quark star of mass M = 0.7M Q with B = 471.3 
MeV fm -3 . Because of the large difference in compactness (or average density) of the quark star, the set of frequencies 
and damping rates of these QNMs differ considerably from that for quark stars with conventional values for B. This 
implies that gravitational wave observations have strong potential to test Nakamura's formation scenario for quark 
stars. In the near future, we can get information about the EOS of quark matter by detecting gravitational waves 
from objects composed of quark matter. Based on the information obtained, we can test the current understanding 
of quark matter and may obtain a more realistic picture of cold quark matter and nonperturbative QCD. 
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TABLE I: Properties of quark stars constructed by the bag model EOS. 





B (MeV/fm 3 ) 


Raa (km) 


Pc (g/cm 3 ) 


M/Mq 


R (km) 


M/R 


A\ 


28.90 


3.800 


2.090 x 


10 14 


2.329 x 


10" 


2 


3.765 


9.135 x 10- 


3 


A2 


28.90 


6.000 


2.135 x 


10 14 


8.877 x 




2 


5.864 


2.235 x 10- 


2 


A3 


28.90 


8.200 


2.205 x 


10 14 


2.165 x 


io- 


1 


7.859 


4.067 x 10- 


2 


Bl 


56.00 


3.800 


4.104 x 


10 14 


4.422 x 


10" 


2 


3.733 


1.749 x 10- 


2 


B2 


56.00 


6.000 


4.285 x 


10 14 


1.637 x 


io- 


1 


5.742 


4.210 x 10- 


2 


B3 


56.00 


8.200 


4.585 x 


10 14 


3.839 x 


io- 


1 


7.560 


7.498 x 10- 


2 


CI 


94.92 


3.800 


7.095 x 


10 14 


7.283 x 


io- 


2 


3.688 


2.916 x 10- 


2 


C2 


94.92 


6.000 


7.662 x 


10 14 


2.592 x 


10" 


1 


5.573 


6.867 x 10- 


2 


C3 


94.92 


8.200 


8.730 x 


10 14 


5.776 x 


io- 


1 


7.156 


1.192 x 10- 


1 


ATI 


471.3 


5.661 


1.486 x 


10 16 


7.000 x 


10" 


1 


3.856 


2.681 x 10- 


1 



TABLE II: Comparison between the /-mode QNMs of quark stars and value given by the empirical formulas 14.111 and 14.211 . In 
this table, Re(u) e m P and /m(w) em p express the range of frequency and damping rate respectively given by using the empirical 
formulas. 





Re(uj) (kHz) 


Im{ui) (Hz) 


Re(uj) em p (kHz) 


Im(uj) emp (Hz) 


Al 


6.822 


2.98 x 10- 


3 


10.52 - 


10.72 


(2.97- 3.42) x 10" 2 


A2 


6.872 


2.67 x 10- 


2 


10.55 - 


10.74 


(2.66-3.10) x 10" 1 


A3 


6.945 


1.12 x 10- 


1 


10.58 - 


10.78 


1.11 - 1.32 


Bl 


9.540 


2.05 x 10- 


2 


12.77- 


12.99 


(2.04 - 2.37) x 10 _1 


B2 


9.676 


1.69 x 10- 


1 


12.84 - 


13.06 


1.67 - 2.00 


B3 


9.873 


6.40 x 10- 


1 


12.95 - 


13.17 


6.17 - 7.74 


CI 


12.50 


9.18 x 10- 


2 


15.20 - 


15.46 


0.913 - 1.07 


C2 


12.80 


6.83 x 10- 


1 


15.36 - 


15.62 


6.63-8.24 


C3 


13.26 


2.28 




15.64 - 


15.90 


20.4 - 28.1 


Nl 


38.10 


20.6 




34.90 - 


35.39 


22.5 - 251.3 
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TABLE III: Deviation from the new empirical formula 14,31 for each stellar model in FigEI I n this table, B emp is the value of 
the bag constant which is calculated by using the empirical formula Ij4,3^ obtained in Sec. IIVI 



Roo (km) 


B (MeV/fm 3 ) 


Re(tu) (kHz) 


B emp (MeV/fm 3 ) 


deviation (%) 


3.8 


28.90 


6.822 


28.67 


-0.813 


3.8 


42.00 


8.242 


41.02 


-2.344 


3.8 


56.00 


9.540 


54.10 


-3.385 


3.8 


75.00 


11.076 


71.79 


-4.285 


3.8 


94.92 


12.503 


90.30 


-4.868 


6.0 


28.90 


6.872 


29.07 


0.581 


6.0 


42.00 


8.331 


41.85 


-0.355 


6.0 


56.00 


9.676 


55.57 


-0.768 


6.0 


75.00 


11.286 


74.39 


-0.819 


6.0 


94.92 


12.803 


94.44 


-0.503 


8.2 


28.90 


6.945 


29.65 


2.608 


8.2 


42.00 


8.458 


43.07 


2.554 


8.2 


56.00 


9.873 


57.74 


3.101 


8.2 


75.00 


11.597 


78.31 


4.415 


8.2 


94.92 


13.257 


100.89 


6.294 
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FIG. 1: Relations (a) between mass M and central density p c , (b) between M and radius R (solid line) and between M and 
radiation radius Roo (broken line) , and (c) between "average density" p = ZM/A-kR 3 and radiation radius R x ol quark stars 
for four values of the bag constant, B = 28.9MeV fm" 3 , 56.0MeV fm^ 3 , 94.92MeV fm~ 3 and 471.3MeV fm~ 3 . 
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FIG. 2: Complex frequencies of /-mode for quark stars. The squares, triangles and circles correspond to stars with B = 28.9 
MeV fm -3 , 56.0 MeV fm -3 , and 94.92 MeV fm -3 , respectively. In each set, the upper, middle and lower marks are for the 
cases of Roa = 8.2 km, 6.0 km, and 3.8 km, respectively. The labels in this figure correspond to the stellar model in Table [I] 
For (b), the plot for B — 471.3 MeV fm -3 is also included as the diamond. 
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FIG. 3: Complex frequencies of both w- and w n -modes (a) for B = 28.9 MeV fm" 3 , (b) for B = 56.0 MeV fm" 3 , (c) for 
B = 94.92 MeV fm" 3 , and (d) for B = 471.3 MeV fm" 3 . For (a)-(c), the upper, middle and lower sequences correspond to 
the quark stars of Roo — 3.8 km, 6.0 km, and 8.2 km, respectively. For (d), the mass of the stellar model is 0.7Mq. Filled and 
non-filled marks denote w- and uin-modes, respectively. These stellar models are listed in Table Q 
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FIG. 4: Complex frequencies of the lowest wn-mode for each stellar model except for the plot for B = 471.3MeV fm 3 , for 
which the second wn-mode is plotted. The labels in this figure correspond to those of the stellar models in Table 
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FIG. 5: The horizontal axis is the bag constant B (MeV/fm -3 ) and the vertical axis is the /-mode frequency Re(u>) . The 
squares, circles and triangles correspond to the stellar models whose radiation radii are fixed 3.8, 6.0 and 8.2 km, respectively. 
The dotted line denotes the new empirical formula 14.31 obtained in Sec. II VI between Re(uj) of /-mode and B. 



